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ABSTRACT 

The  linked  list  and  indexed  list  future  event  sets  are 
investigated  here.  The  interaction  hold  model  and  the 
Jackson  network  model  are  the  underlying  stochastic 
models  considered.  For  the  interaction  hold  model 
and  for  the  (doubly)  linked  list,  we  find,  for  example, 
the  mean  number  of  key  comparisons  performed  in 
order  to  find  a  record’s  insertion  point  into  the  list; 
this  is  useful  when  deciding  whether  to  scan  from  the 
head  or  the  tail  of  the  list.  The  distribution  of  the 
relative  position  of  the  to-be-inserted  record  is  also 
obtained;  for  indexed  lists  this  is  helpful  when  de¬ 
ciding  the  number  of  sub  lists  and  position  (s)  of  the 
middle  pointer(s).  The  Jackson  network  model  has 
a  realistic  event  logic,  but  events  are  restricted  to 
be  exponentially  distributed.  Because  the  stationary 
probabilities  can  be  computed  for  this  model,  it  is 
then  possible  to  evaluate  and  compare  the  (steady- 
state)  performance  of  certain  future  event  sets  (e.g., 
linked  lists  scanned  from  the  head  or  the  tail). 

1  INTRODUCTION 

In  a  discrete-event  simulation,  the  future  event  set 
(FES),  or  also  simulation  calendar,  is  the  data  struc¬ 
ture  that  holds  the  records  of  the  events  scheduled  to 
occur  in  the  future.  Part  of  the  computational  effort 
in  a  simulation  is  devoted  to  the  manipulation  of  the 
FES  as  records  get  removed,  inserted,  updated  and 
looked  at  (e.g.,  for  checking  certain  conditions).  See 
Devroye  (1986),  Manipulation  of  the  FES  must  be 
efficient,  especially  for  large-scale  simulations.  Be¬ 
cause  of  the  various  demands  a  complex  simulation 
puts  on  the  FES,  it  is  not  evident  that  in  such  a  var¬ 
ied  environment  the  more  modern  algorithms  (e.g., 
splay  trees,  binomial  queues,  etc.)  perform  uniformly 
better  than  the  traditional  ones  (e.g.,  indexed  lists, 
heaps,  etc.).  Jones  (1986)  contains  an  extensive  em¬ 
pirical  investigation. 


It  is  clearly  of  interest  to  study  the  performance 
of  the  classical  data  structures  (e.g.,  linked  lists)  in 
a  simulation  environment,  but  the  problem  is  diffi¬ 
cult.  There  are  two  general  approaches  to  analyz¬ 
ing  the  performance  of  a  FES.  The  first  one  is  to 
run  various  real-world  simulations  and  perhaps  arrive 
to  conclusions  based  upon  the  empirical  experiments. 
The  other  approach  is  to  assume  that  the  underlying 
stochastic  system  is  simple  enough  so  that  an  analysis 
of  the  FES  performance  is  tractable.  This  is  the  ap¬ 
proach  taken  here.  The  two  stochastic  models  consid¬ 
ered  are  the  interaction  hold  model  and  the  Jackson 
network  model.  The  former  has  a  simple  event  logic 
but  allows  for  general  probability  functions,  while  the 
latter  has  a  more  realistic  event  logic,  but  events  are 
restricted  to  be  exponentially  distributed.  We  will  be 
looking  at  (doubly)  linked  lists  and  indexed  lists  only. 

Under  the  hold  model  assumption,  first  considered 
for  FES  performance  evaluation  by  Vaucher  (1977), 
the  underlying  stochastic  process  consists  of  a  fixed 
number  L  of  independent  renewal  processes  with 
identical  event  lifetime  distributions.  Each  time  an 
event  triggers  a  transition  (at  the  end  of  its  lifetime), 
its  record  is  removed  from  the  list,  the  event  is  then 
generated  anew  and  its  record  inserted  back  into  the 
FES.  The  hold  model  where  events  may  have  differ¬ 
ent  distribution  functions  was  considered  by  McCor¬ 
mack  and  Sargent  (1981),  and  was  coined  the  interac¬ 
tion  hold  model.  A  multi-server  queue  with  possibly 
nonequivalent  servers  and  with  an  infinite  supply  of 
customers  readily  available  fits  the  interaction  hold 
model  framework.  But,  as  previously  mentioned,  the 
interaction  hold  model  does  have  a  simple  event  logic. 
For  example,  the  FES  will  have  constant  size  L.  For 
the  Jackson  network  model,  introduced  in  Damerdji 
and  Glynn  (1995),  it  is  possible  to  numerically  com¬ 
pute  the  stationary  probabilities,  and  one  can  then 
evaluate  the  performance  of  certain  data  structures 
(e.g.,  linked  list)  under  this  model. 

In  a  simulation,  the  FES  is  typically  kept  sorted. 


based  upon  the  events’  occurrence  times.  Whenever  a 
record  is  to  be  inserted  into  the  FES,  a  search  is  per¬ 
formed  to  find  its  insertion  point  (so  the  FES  stays 
sorted).  The  search  is  carried  out  by  comparing  the 
to-be-inserted  event’s  occurrence  time  with  the  occur¬ 
rence  times  of  the  successive  events  in  the  FES  until 
the  correct  insertion  point  is  found.  The  removal  of 
the  current  event’s  record  and  the  actual  update  of 
the  FES  are  fixed  overheads.  Algorithm  efficiency 
is  measured  here  in  terms  of  number  of  comparisons 
performed  to  find  the  insertion  point  at  a  transition. 

For  the  interaction  hold  model,  the  clock-vector 
across  transitions  is  viewed  as  a  general  state  space 
Markov  chain,  whose  stationary  distribution  can  be 
derived  (see  Section  2).  The  mean  of  the  number 
of  comparisons  per  transition  for  a  linked  list  is  pro¬ 
vided,  also  in  Section  2.  It  is  found,  for  example,  that 
if  the  events  are  all  exponentially  distributed,  it  is 
strictly  better  to  scan  from  the  head  of  the  list  (unless 
the  rates  are  all  equal  in  which  case  it  is  equivalent 
to  scan  from  the  tail  of  the  list).  Indexed  lists,  i.e., 
lists  with  sublists  or  also  buckets,  are  investigated  un¬ 
der  the  interaction  hold  model  assumption.  Here,  we 
only  consider  sublists  that  contain  fixed  numbers  of 
records.  In  the  case  of  an  indexed  list  with  two  sub¬ 
lists  and  with  exponentially  distributed  events,  ex¬ 
plicit  expressions  for  the  number  of  comparisons  per 
transition  in  terms  of  the  event  rates  are  given;  see 
Section  2.  Section  3  contains  an  example. 

The  Jackson  network  model  is  discussed  in  Sec¬ 
tion  4.  An  expression  for  the  expected  number  of 
comparisons  per  transition  in  terms  of  the  stationary 
probabilities  in  given  for  a  general  FES.  An  exam¬ 
ple  illustrates  our  point.  The  Jackson  network  model 
is  introduced  in  the  hope  that  something  might  be 
learned  (e.g.,  is  it  better  to  scan  a  linked  list  from 
the  head  rather  than  the  tail?),  and  then  that  the 
exponential-case  conclusions  might  carry  over  to  the 
general  case,  at  least  to  some  extent.  A  limited  empir¬ 
ical  investigation  is  presented  in  Section  4.  Section  5 
is  the  conclusion. 

2  THE  INTERACTION  HOLD  MODEL 

The  interaction  hold  model  is  now  described  in  more 
detail.  We  will  have  clocks  keep  track  of  the  events’ 
residual  lifetimes,  i.e.,  the  times  remaining  for  the 
events  to  trigger  the  transition.  For  each  event  z,  let 
(7i(.)  be  its  distribution  function  and  (r^j  :  j  >  1) 
its  sequence  of  lifetimes.  Let  =  (^n,i)  •  •  • 
denote  the  clock-vector  reading  at  the  n’th  transition 
(in  fact,  at  the  instant  after  the  transition).  The  event 
that  will  trigger  the  next  transition  has  the  smallest 
clock  reading.  If  event  z,  say,  triggers  the  n’th  transi¬ 


tion,  the  clocks  corresponding  to  the  other  events  all 
get  updated,  while  event  z  gets  generated  anew  from 
G'i(.).  Its  clock  reading  will  be  where  Ni{n) 

is  the  number  of  lifetimes  of  event  z  that  have  been 
observed  by  the  n’th  transition.  We  have  that 

L 

Xfi—  ^  ^  I  n— l,i  ~  j  ^-^n-1,1 

t=l  ^  ’  *’ 

— An-l,i,  .  .  *  ,  Ti^Ni{n)^  •  •  •  j  » 

where  /[.]  is  the  indicator  event.  Since  the  triggering 
event  is  generated  anew  independently  of  the  past,  it 
is  clear  that  {Xn  :  n  >  0}  is  a  Markov  chain,  with 
state  space  =  {(vi,...,yL)  :  yi  >  0, 1  <  z  < 
L}.  It  is  shown  in  Damerdji  and  Glynn  (1995)  that 
this  general  state  space  Markov  chain  is  recurrent  in 
the  sense  of  Harris  (see  Meyn  and  Tweedie  (1993)  for 
a  definition),  and  that  its  stationary  distribution  is 
given  by 


7r(dxi, 


L 

■  ■ .  d^L)  = 

t=i 


Ai 


Gi{dxi) 


^3  ^3  (^3  )  ’ 

j^i 


where  Gj{.)  =  1  —  Gj{.).  Let  be  the  probability 
measure  of  the  Markov  chain  on  its  infinite  path  space 
with  initial  distribution  tt. 

As  previously  mentioned,  efficiency  of  a  FES  is 
measured  here  in  terms  of  number  of  comparisons  (in 
order  to  find  the  correct  insertion  point)  per  scan. 
Doubly  linked  lists  are  considered  first.  Let  be 
the  random  variable  that  denotes  the  relative  posi¬ 
tion  (minus  one)  of  the  event-record  inserted  at  the 
/c’th  transition.  Let  also  (respectively,  Fb,A;)  rep¬ 
resent  the  number  of  comparisons  (minus  one)  per¬ 
formed  at  the  A:’th  transition  when  scanning  from  the 
head  (resp.,  tail)  of  the  list.  For  example,  if  at  the 
A;’th  transition,  the  generated  event  has  the  second 
smallest  time,  then  Wk  =  1,  Ff,*;  =  1,  and  = 
L  —  1  —  =  L  —2.  We  obtain  the  following. 

Proposition  1, 


Gi{t)Gj{t)dt, 


where  is  the  expectation  under 

This  result  is  a  simplified  version  of  an  equation 
that  appears  in  McCormack  and  Sargent  (1981).  For 
identical  distributions,  see  Vaucher  (1977).  Our  goal 
now  is  to  investigate  whether  it  is  more  efficient  to 


scan  the  list  from  the  head  or  the  tail.  Certain  prop¬ 
erties  of  the  distributions  may  be  determining  in  some 
cases.  This  is  discussed  next. 

Definition  1,  A  distribution  function  G,  with  finite 
mean  1/A,  is  said  to  be  new  better  (resp.f  worse)  than 
used  in  expectation  (NBUE  (NWUE))  if 

G{t)dt  <  (>)  (l/A)G(a:)  for  all  a:  >  0. 


FVom  Barlow  and  Proschan  (1975),  it  follows  that 
if  two  distribution  functions  Gi  and  Gj,  with  means 
1/Ai  and  1/Aj,  are  NBUE  (resp.,  NWUE),  then 

rGiit)cjit)dt  >  (<) 

Jo 

Because  the  expression  in  Proposition  1  contains  such 
an  integral,  the  following  corollary  ensues. 

Corollary  1.  If  the  event  lifetimes  are  all  NBUE 
(resp.y  NWUE)f  then 


E^Yra  > 


(<) 


1 


I 


EE 

1=1 


XiXj 

Xi  d"  Xj 


When  the  distributions  are  all  NBUE,  the  lower 
bound  of  Corollary  1  complements  McCormack  and 
Sargent  (1981),  who  provide  an  upper  bound  (for  gen¬ 
eral  distributions).  Consider  the  following  lemma. 

Lemma  1,  For  any  positive  and  finite  numbers  Ai, 
. . .  ,  Ax,,  we  have  that 


1 


L 


EE 

z=l  ji^i 


X{Xj 
At  -h  Aj 


< 


L  -  1 
2 


with  equality  if  and  only  if  Xi  =  . . .  =  A/,. 

An  immediate  consequence  of  Proposition  1  and 
Lemma  1  is  the  following. 

Corollary  2.  If  the  events  are  all  exponentially  dis- 
tributedj  it  is  strictly  better  to  scan  from  the  head  of 
the  listj  unless  the  events  have  identical  rates  in  which 
case  it  is  equivalent  to  scan  from  the  back. 

Because  the  stationary  distribution  tt  of  the 
Markov  chain  is  known,  it  is  possible  to  compute  the 
distribution  =  •]  of  the  relative  position  of 

the  record  to  be  inserted.  See  Damerdji  and  Glynn 
(1995)  for  a  general  expression.  If  the  events  are  all 
exponentially  distributed,  the  distribution  of  de¬ 
pends  only  upon  the  event  rates.  Let  C{g,h)  be  the 
number  of  possible  combinations  of  size  h  out  of  a 
set  of  size  p,  Qi{j,n)  be  the  fth  set  of  size  n  out 


of  (1,...,L}  “  {i},  Rj{i,n)  be  its  complement  in 
{1, . . . ,  L}  -  {i},  and  Tb,a,q{iy  j,  n)  be  the  6’th  subset 
of  size  a  out  of  —  {^}. 

Proposition  2.  If  the  events  are  all  exponentially 
distributed^  we  have,  for  n  =  0, . . . ,  L  —  1,  that 


L  .2  <7(1/— l,n) 

=  n]  =  ^  X) 

i=i  *  i=i 

veRj{i,n)  9€Qj(i,n)  o=0 

C(n— l,a) 

•[  E  E 

6=1 

+  E  . 

Also,  we  have  that  P»[Wi  =  0]  = 


In  an  indexed  list,  i.e.,  a  list  with  sublists,  the 
middle  pointers  point  to  certain  records  (possibly 
dummy  records)  in  the  list  chosen  depending  upon 
the  strategy  followed.  One  approach  is  to  keep  the 
sublists  with  constant  sizes.  Another  strategy  is  to 
have  the  middle  records  chosen  based  upon  their 
event  times  of  occurrence.  This  was  the  setting  of 
Engelbrecht-Wiggans  and  Maxwell  (1978)  and  Davey 
and  Vaucher  (1980).  Here,  the  former  approach  is  in¬ 
vestigated,  and  so  each  sublist  will  have  a  fixed  num¬ 
ber  of  records. 

For  simplicity  of  analysis,  consider  now  an  indexed 
list  with  two  sublists  only.  Let  m  be  the  middle 
record,  i.e.,  the  record  such  that  the  (unique  here) 
middle  pointer  is  pointing  to.  Because  of  the  two- 
way  links,  four  search  strategies  are  possible:  the  first 
(resp.  second)  sublist  can  be  scanned  forward  from 
the  head  (resp.  middle)  pointer  or  backward  from  the 
middle  (resp.  tail)  pointer.  Let  ZFF,*;(7n), 

^BF,fc(^)j  and  ZBB,k(^)  be  the  number  of  compar¬ 
isons  (minus  one)  performed  at  the  fc’th  transition  in 
order  to  find  the  correct  insertion  point  for  the  differ¬ 
ent  strategies.  Explicit  expressions  can  be  obtained. 
For  example. 


m— 2  L—1 

=  1+  ^  na(n)  +  E  (n  — m-l-l)a(n), 

n=0  n=m— 1 

where  a(n)  =  =  w]-  See  Damerdji  and  Glynn 

(1995)  for  the  other  three  strategies.  It  is  then  possi¬ 
ble  to  optimize  with  respect  to  m,  and  also  compare 


the  performances  for  the  four  strategies.  An  example 
is  given  in  the  next  section.  Lists  with  more  than  two 
sublists  can  be  similarly  investigated. 

3  AN  EXAMPLE 

We  consider  here  an  interaction  hold  model  with  ten 
exponentially  distributed  events,  split  into  two  equal- 
size  groups:  events  1  thru  5  have  rates  Ai  =  ...  = 
A5  —  1.0,  while  events  6  thru  10  have  rates  Xq—  . . .  = 
A 10.  Three  experiments  are  conducted:  Ae  is  chosen 
such  that  l/Ae  =  LO  in  the  first  one  (the  two  groups 
are  equal),  l/Ae  —  2.0  in  the  second  one  (the  two 
groups  are  rather  different),  and  l/Ae  =  10.0  in  the 
last  experiment  (the  two  groups  are  veiy  different). 


l/Ae 

1.0 

2.0 

10.0 

4.500 

4.222 

2.826 

4.500 

4.778 

6.174 

Table  1:  Linked  list. 


For  the  third  experiment,  i.e.,  when  the  two  groups 
are  very  different,  it  is  seen  from  Table  (1)  that  it  is 
far  better  to  scan  the  list  from  the  head  pointer.  From 
Table  (2),  the  probabilities  of  the  relative  position  of 
the  inserted  record  are  equal  in  the  first  experiment 
because  of  the  memoryless  property  of  the  exponen¬ 
tial  distribution  and  by  symmetiy.  Note  that  the 
probabilities  are  decreasing  (relative  to  the  position) 
for  the  second  and  third  experiments.  For  the  third 
experiment  for  example,  it  is  almost  ten  times  more 
likely  that  the  newly  inserted  record  be  at  the  top 
of  the  list  rather  than  the  bottom  (0.166  vs.  0.018). 
This  stresses  again  that,  at  least  in  the  exponential 
case,  it  is  better  to  scan  a  linked  list  from  the  head 
rather  than  the  tail. 

Once  the  distribution  of  Wi  is  obtained,  it  is  strai¬ 
ghtforward  to  compute 

^7r^BF,i(^),  and  Er^ZBE^iim)  for  m  =  3, . . . ,  9,  and 
optimize  over  m  for  each  strategy.  Let  mpp,  mpg, 
TTigp,  and  mgg  be  the  respective  optimum.  For 
the  second  and  third  experiments,  it  is  found  that 
the  best  strategy  is  forward- forward,  closely  followed 
by  backward-forward.  Backward-backward  was  the 
worst  strategy.  For  the  third  experiment  and  for  the 
forward- forward  strategy,  the  fourth  record  is  the  op¬ 
timal  middle  record;  in  that  case  the  first  sublist  con¬ 
tains  three  records,  while  the  second  sublist  contains 
seven  records.  The  first  sublist  has  size  smaller  than 
50%  that  of  the  second  one.  Comparing  the  best 
performance  for  indexed  list  and  regular  linked  list 
for  the  third  experiment  (E7r^FF,i(^)  =  2.299  and 


l/Ae 

1.0 

2.0 

10.0 

a(0) 

.100 

.111 

.166 

a(l) 

.100 

.109 

.164 

a(2) 

.100 

.108 

.160 

a(3) 

.100 

.106 

.153 

a  (4) 

.100 

.104 

.140 

a(5) 

.100 

,101 

.095 

a{6) 

.100 

.098 

.052 

a(7) 

.100 

.093 

.028 

a(8) 

.100 

.087 

.020 

a(9) 

.100 

.078 

.018 

Table  2:  Distribution  ofWi. 


l/Ae 

1.0 

2.0 

10.0 

TTL  pp 

6 

6 

4 

ETrZFF,lini*) 

3.000 

2.924 

2.299 

mpB 

6 

6 

6 

E-^EFB^im*) 

3.000 

3.040 

3.120 

rngp 

6 

6 

4 

ETrZBF,l{m*) 

3.000 

2.957 

2.312 

6 

6 

6 

3.000 

3.073 

3.248 

Table  3:  Indexed  list  with  two  sublists. 


=  2.826),  it  is  not  evident  that  the  indexed 
list  far  outperforms  the  regular  list  given  the  overhead 
associated  with  maintaining  the  indexed  list. 

4  THE  JACKSON  NETWORK  MODEL 

Open  and  closed  Jackson  queueing  networks  are  pro¬ 
posed  as  stochastic  models  for  FES  performance  eval¬ 
uation.  As  previously  mentioned,  these  queueing  net¬ 
work  systems  have  a  reasonably  complex  event  logic, 
although  only  exponentially  distributed  event  times 
are  considered.  When  simulating  a  Markovian  sys¬ 
tem  such  as  a  Jackson  network,  one  need  not  even 
use  a  calendar  (from  the  properties  of  the  exponen¬ 
tial  distribution).  This  class  of  stochastic  models  is 
proposed  in  the  hope  that  what  might  be  learned  in 
the  exponential  case  extends  to  the  general  case,  at 
least  to  some  extent. 

Only  open  networks  are  discussed  here.  There  are 
d  single-server  stations  with  external  arrival  rate  A^  to 
station  i.  Service  at  station  i  is  with  rate  fxi  and  is  on 
a  first-in-first-out  basis.  Upon  service  completion,  the 
customer  departs  the  system  with  probability  pio  or 


is  routed  to  station  j  with  probability  pij.  The  state 
of  the  system  can  be  described  by  u  =  (ui, . . .  .Ud), 
where  Ui  is  the  number  of  customers  at  station  i.  If 
Ck  is  the  A:-unit  vector,  the  state  changes  from  u  to 
tl  -f  Ci  if  an  external  arrival  to  station  i  occurs,  from 
u  to  u  —  e*  if  a  departure  from  station  i  to  the  outside 
system  occurs,  and  from  n  to  iT  —  e,  +  if  an  end  of 
service  at  station  i  occurs  and  the  departing  customer 
is  routed  to  station  j. 

Let  x(-)  the  stationary  distribution,  which  is 
given  by 

d 

x(«)  =  JJ(1  -  Oiilni){ai/ni)''‘, 

t=i 

where  a  =  (ai, . . . ,  q^)  is  a  vector  solution  to  the  traf¬ 
fic  equations.  Let  Ai  and  Si  be  a  station  i  generic  ex¬ 
ternal  interarrival  time  and  service  time,  respectively. 
For  any  state  w,  let  b{w)  be  the  set  of  busy  servers  in 
that  state,  A  general  expression  for  the  steady-state 
expected  number  of  comparisons  per  transition  for  a 
general  FES  can  be  provided  by  conditioning  on  the 
first  two  states  of  the  system.  We  have  that 

E[#  comparisons/transition]  = 

w  £=1  i^b{w) 

.  s  j]  f^iPijXi'^) 

u:ui>2;uj>l  v=u  —  ei+ej 

+  E  E  E[#|5i,  then  Sjlmpijxiu) 

•u:ui>2;uj=0  v—u  —  ei+Cj 

+  E  E 

u:Ui>2  v—u  —  Ci 

+  E  E  E[#|Aj,then  Sj]Aix(w) 

n:ui=0  v=u-hei 

+  E  E  E[#M<]Aix(tr)}. 

v=u-\-ei 

We  are  assuming  here  that  if  an  external  arrival  oc¬ 
curs  to  an  empty  station,  the  next  interarrival  time  is 
generated  and  inserted  into  the  FES  before  the  service 
time.  Also,  when  a  departure  from  a  station  with  a 
nonempty  queue  to  an  empty  station  occurs,  the  ser¬ 
vice  time  of  the  first  awaiting  customer  is  generated 
(and  inserted)  before  the  service  time  of  the  arriving 
customer. 

As  previously  mentioned,  it  is  possible  to  compute 
the  stationary  probabilities.  In  order  to  evaluate  the 
performance  of  a  FES,  it  suffices  then  to  compute  the 


A 

1.0 

1.0 

1.0 

1.0 

1.0 

Pi 

1.2 

2.0 

10.0 

2.0 

10.0 

P2 

1.2 

2.0 

10.0 

3.0 

5.0 

F 

1.883 

1.611 

1.154 

1.522 

1.216 

B 

1.893 

1.722 

1.645 

1.700 

1.650 

Eio% 

1.982 

1.498 

1.000 

1.000 

1.000 

Bio% 

1.350 

1.500 

1.666 

1.666 

1.666 

■^50% 

2.182 

1.490 

1.000 

1.319 

1.000 

■650% 

1.399 

1.527 

1.666 

1.574 

1.666 

Table  4:  Linked  list,  open  Jackson  network. 


inner  expected  values.  A  similar  analysis  for  closed 
networks  follows  (see  Damerdji  and  Glynn  (1995)). 
An  example  for  an  open  network  and  for  the  linked 
list  is  now  given. 

Consider  a  two-station  in  tandem  with  external  ar¬ 
rivals  to  the  first  station  only.  It  follows  then  that 
A2  =  0  and  Pi2  =  P20  ~  1-  Let  A  be  the  arrival  rate 
to  station  1,  fixed  to  1.  For  the  stationary  probabili¬ 
ties  to  exist,  it  must  be  that  pi  >  I  and  /i2  >  L  For 
all  the  values  of  pi  and  p2  tried,  we  find  that  it  is 
better  to  scan  from  the  head  of  the  list  rather  its  tail. 
See  Table  (4),  where  F  (resp.  B)  denotes  the  expected 
number  of  comparisons  per  transition  when  scanning 
from  the  head  (resp.  tail)  of  the  list.  For  example, 
for  /ij  =  /i2  =  10,  it  is  far  better  to  scan  from  the 
head  rather  than  the  tail  {F  ~  1.154  and  B  —  1.645). 
When  the  rates  are  closer  to  one  another,  it  is  only 
slightly  better  to  scan  from  the  head. 

We  now  empirically  investigate  whether  results  for 
the  exponential  case  do  extend  to  the  nonexponential 
one.  Uniform  distributions  are  considered  in  the  two 
experiments  performed.  If  r  is  the  rate  of  an  event 
in  the  exponential  case  (and  so  1/r  is  the  mean), 
the  uniform  distributions  U{a,b)  for  that  event  are 
such  b  —  a  —  10%  of  the  mean  in  the  first  experiment 
and  50%  of  the  mean  in  the  second.  The  tandem- 
queue  was  then  simulated  over  one  million  transi¬ 
tions.  Let  Fio%,  Biq%,  ^50%,  and  ^50%  be  the  re¬ 
spective  estimates  of  the  expected  number  of  compar¬ 
isons  per  scan.  See  Table  (4)  for  the  results.  When 
the  rates  are  all  close  to  one  another  (e.g.,  A  =  1, 

=  ;x2  =  1.2),  it  is  better  to  scan  from  the  back  of 
the  list.  We  note  that  when  it  is  far  more  efficient  to 
scan  from  the  head  in  the  exponential  case,  it  is  also 
better  to  scan  from  the  head  in  the  nonexponential 
case  as  well.  When  it  is  not  far  more  efficient  to  scan 
from  the  head  in  the  exponential  case,  we  note  that 
it  is  better  to  scan  from  the  back  of  the  list  in  the 
nonexponential  case. 


5  CONCLUSION 

Performance  evaluation  of  linked  lists  and  indexed 
lists  has  been  investigated  here.  The  stochastic  mod¬ 
els  considered,  namely  the  interaction  hold  model  and 
the  Jackson  network  model  are  very  different.  The 
former  has  a  simple  event  logic  but  allows  for  general 
distribution  functions,  while  the  latter  has  a  more 
complex  event  logic,  but  distribution  functions  are 
exponential.  Our  analytic  and  empirical  results  sug¬ 
gest  that  in  a  simulation  with  different  types  of  dis¬ 
tribution  functions  with  the  short  events  frequently 
regenerated  upon  transition,  it  is  often  better  to  scan 
a  linked  list  from  the  head  rather  than  from  the  tail. 
However,  there  will  also  be  situations  where  it  is  bet¬ 
ter  to  scan  from  the  back  of  the  list.  For  example, 
if  the  distributions  are  Uniform  (100,120),  IViangu- 
lar  (90,95,110),  and  Normal  (120,5),  it  would  seem 
better  to  scan  from  the  back  of  the  list. 
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